Metastability, Nucleation, and Noise-Enhanced Stabilization Out of Equilibrium 



O 
O 
(N 

> 

O 



o 

a 

CO 



T3 

O 
> 

in 
in 

^J- 

o 

^D 
O 

-i— > 

c3 



i 

c 
o 
o 



X 



Pablo I. Hurtado, 1 J. Marro, 2 and P.L. Garrido 2 

1 Laboratoire des Colloides, Verves et Nanomateriaux, 

Universite Montpellier II, Montpellier 34095, France 

^Institute Carlos I for Theoretical and Computational Physics, 

and Departamento de Electromagnetismo y Fisica de la Materia, 

Universidad de Granada, E- 18071 -Granada, Espana 

We study metastability and nucleation in a kinetic 2D-Ising model which is driven out of equi- 
librium by a small random perturbation of the usual dynamics at temperature T. We show that, at 
a mesoscopic/cluster level, a nonequilibrium potential describes in a simple way metastable states 
and their decay. We thus predict noise-enhanced stability of the metastable phase and resonant 
propagation of domain walls at low-T. This follows from the nonlinear interplay between thermal 
and nonequilibrium fluctuations, which induces reentrant behavior of the surface tension as a func- 
tion of T. Our results, which are confirmed by Monte Carlo simulations, can be also understood in 
terms of a Langevin equation with competing additive and multiplicative noises. 
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Relaxation in many natural systems proceeds through 
metastable states 0, 0, i, H |E IS S II II- This is of- 
ten observed in condensed matter physics itg|, and also 
in various other fields, from cosmology 6] to biology 
[7J and high-energy physics [8j. In spite of such ubiq- 
uity, the microscopic understanding of metastability still 
raises fundamental questions. A main difficulty is that 
this is a dynamic phenomenon not included in the en- 
semble formalism. Even so, metastable states may be 
regarded in many cases as equilibrium states for times 
short compared with their relaxation time, and one may 
derive macroscopic properties from restricted ensembles 
[l(, or obtain fluctuation-dissipation theorems [9|. On 
the long run, however, metastable states eventually de- 
cay triggered by rare fluctuations. This relaxation can 
be described in terms of free-energy differences as far 
as one is dealing with systems evolving towards equilib- 
rium steady states [3|. However, as a rule, natural sys- 
tems are open to the environment, which induces currents 
of macroscopic observables or competitions between op- 
posing tendencies which typically break detailed balance 
[10] . Consequently, in many cases, stationary states are 
not equilibrium states, but are strongly influenced by dy- 
namics, which adds further challenge to the microscopic 
understanding of metastability. 

In this paper we report on the nature of metastabil- 
ity in a full nonequilibrium setting. In order to focus 
on the basic physics, we study the simplest nonequi- 
librium model in which metastable states are relevant, 
namely an Ising model with dynamic impurities. We 
show that introducing a nonequilibrium condition has 
dramatical effects on metastable dynamics. In partic- 
ular, we find noise-enhanced stabilization of metastable 
states, resonant propagation of domain walls and other 
novel low temperature physics not observed in equilib- 
rium. Surprisingly, these nonequilibrium phenomena can 
be understood via extended nucleation theory, starting 



from a nonequilibrium potential or free-energy at a meso- 
scopic/cluster level. This is possible because the excita- 
tions responsible for the metastable-stable transition re- 
sult from the competition between a surface and a bulk 
term, as in equilibrium. We thus anticipate that a sim- 
ilar approach can be used to understand metastability 
in many other nonequilibrium systems in which a sur- 
face/bulk competition is a dominant mechanism. Our 
results are also relevant for metastable nanodevices where 
nonequilibrium impurities play a fundamental role [10j. 
Consider a two-dimensional square lattice of side L 
with periodic boundary conditions. We define a spin vari- 
able Si = ±1 at each node, i G [1, N = L 2 ]. Spins interact 
among them and with an external magnetic field h via the 
Ising Hamiltonian function, H = — J2uj) s i s j ~ hJ2i s *i 
where the first sum runs over all nearest-neighbors pairs. 
We also define a stochastic single spin-flip dynamics with 
transition rate 



w(s -► s 4 ) = p + (1 - p)y(pAHi), 



(1) 



where s = {sk} and s l stand for the configurations before 
and after flipping the spin at node i, respectively, A7ij 
is the energy increment in such flip, and [3 = 1/T. The 
function in (JTJ) is chosen here as ^(x) = (l + e x )^ 1 , which 
corresponds to the Glauber rate. However, similar results 
hold also for the Metropolis rate, and possibly for many 
other local, detailed-balanced dynamical rules \1/. 

For any < p < 1 two different heat baths com- 
pete in ([1]): One is at temperature T, which operates 
with probability (1 — p), and the other induces com- 
pletely random spin-flips (as a bath at infinite tem- 
perature would do) with probability p. As a result of 
this competition, a nonequilibrium steady state sets in 
asymptotically [10|, [ll|, [l2|, |l3[, which cannot be char- 
acterized by any Gibbsian measure. For h = 0, the 
model exhibits an order-disorder continuous phase tran- 
sition at temperature T c (p) < T c {p = 0) = T ons , where 



Tons = 2/ln(l + \/2), and the order washes out for any 
p > p c = (y/2 - l) 2 w 0.17, even at T = 0. This phase 
transition belongs to the Ising universality class, a result 
that has lead to the belief that two-temperature nonequi- 
librium Ising models behave in general as their equilib- 
rium (p = 0) counterpart. We show below that, in what 
concerns metastability, essential differences indeed exist. 

We remark that our motivation here for the random 
perturbation p is that it is the simplest microscopic mech- 
anism that induces nonequilibrium behavior. However, 
similar dynamic impurities play a fundamental role in 
some natural systems. As a matter of fact, an equiva- 
lent mechanism has been used to model the macroscopic 
consequences of rapidly-diffusing local defects [lfj and 
quantum tunneling [111 ] in magnetic samples and, more 
generally, the origin of scale-invariant avalanches [l4j . 

For small h < and T < T c (p), an initial homoge- 
neous state with all spins up is metastable. It will even- 
tually decay toward the stable state of magnetization 
m = TV -1 Y2i s i < 0- Inspection of escape configura- 
tions (Fig. [T]) shows that this is a highly inhomogeneous 
process triggered by (large) compact clusters of the sta- 
ble phase. These excitations then grow or shrink in the 
metastable sea depending on the competition between 
their surface, which hampers cluster growth, and their 
bulk, which favors it. In equilibrium (p = 0), this compe- 
tition is controlled by the cluster interfacial free-energy, 
or surface tension. Far from equilibrium (0 < p < 1), de- 
spite the lack of a proper bulk free-energy function, one 
may define [13[ an effective surface tension ao(T,p), that 
captures the properties of the nonequilibrium interface. 
This is based on the assumption that the normalization of 
a probability measure for interface configurations can be 
interpreted as a sort of nonequilibrium partition function 





FIG. 1: (Color online) Escape configurations for L = 128, p = 
0.01, h = -0.25 and T/T onB = 0.1, 0.3 and 0.7, respectively 
from left to right. 



FIG. 2: (Color online) Cluster distribution P(C) for T = 
0.2T O ns, L — 53, h = —0.1 and, from bottom to top, p = 0.001, 
0.005, 0.01, 0.02, with 7 = 0.815, 0.82, 0.83, 0.85, respectively. 
Lines are theoretical predictions and points are MC results. 
Inset: 1Z C vs. T ons /T for L = 53, h = —0.1 and, from top to 
bottom, p = 0, 0.001, 0.005 and 0.01. The n th curve (from 
bottom to top) has been rescaled by a factor lO 1 -™ -1 ' (main 
plot), or shifted (4 — n) units (inset) in the y-axis. 



[13J ] . Similar assumptions have been shown to yield excel- 
lent results for nonequilibrium phase transitions [1 51 ] . In- 
teresting enough, o~o(T,p) exhibits non-monotonous tem- 
perature dependence for any < p < 1 with a maximum 
at a non-trivial value of T [l3|, [l6[ . 

Consequently with the surface/bulk competition driv- 
ing the relevant excitations, it is straightforward to write 
down an ansatz for an effective free- energy cost of a clus- 
ter of radius R, T{R) = "f[2nRa -nR 2 2m a \h\}; see @,||]. 
We may then derive the zero-field spontaneous magneti- 
zation mo(T, p) within a mean-field approximation [121 ] 
and the cluster form factor fl(T,p) from Co via the Wulff 
construction [131] : 7 is a weakly-varying parameter, very 
close to one, that stands for small corrections to clas- 
sical nucleation theory [J]. The critical cluster radius, 
such that supercritical (subcritical) clusters tend to grow 
(shrink), thus follows as 1Z C = <Jo/(2mo\h\). Estimates of 
1Z C in MC simulations (see [17| for details) are shown in 
the inset to Fig. [2] together with the analytical predic- 
tions. The agreement is excellent for temperatures well 
below T c (p) and, most important, lZ c (T,p,h) exhibits 
non-monotonous T-dependence for any p > 0. 

Our ansatz above also implies a Boltzmann-like dis- 
tribution, P(C) = M- 1 exp[-/3J r (yc7^)] for the frac- 
tion of stable-phase clusters of volume C — VtR 2 in the 
metastable phase. The normalization M. = 29/(1 — to), 
with 9 = £g=i Cexp[-/3^( v / C7n)] and C* = Wl 2 c , is 
defined so that the metastable state has the mean-field 
magnetization m(T,p,h) [12j. Fig. [2] depicts our results 
for P(C). Again, theoretical predictions compare very 
well to MC results. For p = 0.001, MC data reveal a 
non-trivial structure in P{C) which is not captured by 
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FIG. 3: (Color online) v vs. T ons /T for h = -0.1 and p = 
0, 0.001, 0.005, 0.01, 0.02, 0.03, 0.04, 0.05 (from bottom to 
top). Curves are linear fits to data. Inset: Slope of linear fits 
vs. p. The line is the theoretical prediction, with F « 0.0077. 



our continuous theory. This is due to the lattice struc- 
ture which, for low-T and small p, gives rise to resonances 
for clusters with complete shells, i.e. 2 x 2, 3 x 2, 3 x 3; 
see Fig. [2] For larger p and/or T, fluctuations wash out 
this effect. Also interesting in Fig. is that the nonequi- 
librium perturbation p enhances fluctuations and favors 
larger clusters. 

The assumed F{R) also involves a force per unit area 
which controls the growth of supercritical droplets. In 
particular, the propagation velocity of a large cluster 
should obey the Allen-Calm expression vq = 2r'mo|/i|/<7o 
[2j, where r" is a non-universal constant. Our estimates 
for vo in MC simulations of a flat propagating inter- 
face (an infinitely-large cluster) [l7| are in Fig. [3) The 
fact that ^o exhibits the expected non-monotonous In- 
dependence means that cooling the system favors domain 
wall propagation in a noncquilibrium setting. Moreover, 
an « a(p)T at low-T, with a(p) = ln[(l — y/p)/(jp + y/p)] 
[ll H (l3> so that we expect w « [T/a(p)](T ons /T) 
in this regime, where T is a different non-universal con- 
stant and we assumed too(T — * 0,p) ~ 1. This is nicely 
confirmed in simulations, see inset to Fig. [3l 

The nucleation rate X for critical clusters determines 
the metastable-state lifetime. From our hypothesis, X = 
A\h\ 5 exp[— (3^(71 c )], where A(p) is a non-universal am- 
plitude and S m 3 for random updatings [3]. The re- 
laxation pattern depends on the balance between two 
different length scales, namely, L and the mean cluster 
separation lZo(T,p,h) — (vo/X) 1 / 3 (typically, 1Z C <C 7?.o, 
L). For IZo ^5> L [Single-Droplet (SD) regime], nucle- 
ation of a single critical cluster is the relevant excitation, 
and the metastable-state lifetime is tsd = (T 2 T) _1 . For 
Ko < i [Multidroplet (MD) regime], the metastable- 
stable transition proceeds via the nucleation of many 
critical clusters, and t M d = [^27(3 In 2)] - 1 / 3 '3]. The 
crossover corresponds to the dynamic spinodal point, 
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FIG. 4: (Color online) Lifetime r vs. T for the same p 
and conditions than for the inset in Fig. [2] The n th curve 
(from top to bottom) is rescaled by a factor 10 -2 '™ -1 '. Solid 
lines are theoretical predictions for (from top to bottom) 
7=1, 0.85, 0.77, 0.65. Amplitudes A(p) € [10- 3 ,10~ 2 ]. 
Top inset: |/idsp| vs. T for (from top to bottom) p = 
0, 0.001, 0.005, 0.01, 0.02. Bottom inset: r vs. D, as derived 



from Langevin equation, for h 
increasing from top to bottom. 
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\hT,sp\(T,p) = ^7crg/(6m Tlni), such that the SD 
(MD) regime holds for \h\ < \h BSP \ (\h\ > \h BS p\)- 

We measured t in MC simulations using rejection-free 
methods 17|,[l8|]- This is in Fig. [4] together with our pre- 
dictions. Interesting enough, we observe that r increases 
with T for fixed p > at low-T. That is, the local 
stability of nonequilibrium metastable states is enhanced 
by the addition of thermal noise. This behavior resem- 
bles the noise enhanced stability (NES) phenomenon re- 
ported in experiments on unstable systems [191 ] . which 
is in contrast with the simple Arrhenius law observed 
in equilibrium. Remarkably, only thermal NES is ob- 
served: Increasing p for fixed T always results in shorter 
t. This complex phenomenology is captured by our sim- 
ple ansatz, which traces back the NES phenomenon to 
the low-T anomaly in ero- More generally, the stochas- 
tic resonance in r (and vq) stems from a nonlinear co- 
operative interplay between thermal and nonequilibrium 
noises. Although both noise sources induce disorder 
when applied independently, their combined effect results 
in a resonant stabilization of the metastable phase at 
low-T. This non-linear effect is also reflected in the mor- 
phology of the metastable-stable transition. In particu- 
lar, |/idsp| inherits the non-monotonous T-dependence of 
(TO; see top inset to Fig. [H resulting in a novel MD regime 
at low-T not observed in equilibrium. This is confirmed 
by direct inspection of escape configurations, e.g. Fig. Q] 

One may get further physical insight by rewriting 
the rate (JTJ) as p + (1 - p)^(x t ) = \&(xf), with Xi = 
[3AH.i and xf — /3 c fA7Yi. The resulting effective tem- 
perature T c f(xi,p) = P~ { then follows as T c f/T = 



x % (x i +ln[(l-p)/(l+pe x *)]) \ For any p > 0, T cf 
changes from spin to spin and depends on the local or- 
der (e.g. number of broken bonds for a given spin): the 
smaller the number of broken bonds, the higher the local 
order, and the larger the effective temperature. There- 
fore, for p > 0, the strength of fluctuations affecting a 
spin increases with the local order parameter. This is 
the fingerprint of multiplicative noise, and it allows us 
to write a Langevin equation which captures the essen- 
tial physics. In its simplest, 0-dimcnsional form, this 
equation is dtip = ip — ^ + h + y/D + flip 2 £(i), where 
£(£) is a Gaussian white noise with (£(£)) = and 
(£(*)£(*0) = 2S ( f - *')> D is the strength of the ther- 
mal noise, ft. is a magnetic field, and \x is the renormal- 
ized version of the noncquilibrium parameter p. This 
equation describes a Brownian particle in an asymmet- 
ric bimodal potential, V(tf>) = — \ip 2 + jtp 4 — hip, sub- 
ject to fluctuations which increase with ip 2 and whose 
amplitude remains non-zero as D —> for any \x > 



[13]. A full description of the problem, including 



the compact excitations observed in simulations, would 
of course involve the spatially-extended version of this 
equation. However, the above toy mean-field equa- 
tion already contains the essential competition between 
thermal (D) and nonequilibrium (//) fluctuations in a 
metastable potential which characterizes our system. 
The steady distribution of the stochastically-equivalent 
Fokker-Planck equation in Stratonovich sense is P s t (ip) = 
A[2y/v(D + juffi]- 1 exp[-D- 1 W(i>)], where W(tp) = 
f [ip 2 - (l + rf)ln(rf + ^ 2 )]-/id 1 / 2 tan- 1 (V'^ 1/2 ),with 
d = D/fi and A a normalization constant. The ex- 
trema of the effective potential W(ip) are the same 
than those of V(ip), namely tpk = 2cos(#fc)/v / 3, with 
6 k = |[cos _1 (-|\/t7ft) + 2kn], and k = 0,1,2. For 
h < 0, ipo, ifii and ip2 correspond to the metastable, 
stable and unstable extrema, respectively and the es- 
cape time from the metastable minimum is t{D, \x, h) « 
27:[\V"^o)V"(ijj 2 )\]-^ 2 eMi(W(iP 2 ) - W(ifo))], where 
V" = d 2 V(x). Identifying D with temperature, this ap- 
proach recovers the thermal NES phenomenon for r ob- 
served in the microscopic model (see bottom inset to Fig. 

n. 

Summing up, we have shown that introducing a 
nonequilibrium condition has important effects on 
metastable dynamics. The non-linear interplay between 
thermal and nonequilibrium fluctuations, as captured by 
Co j results in resonant low-T phenomena, e.g., noise- 
enhanced stabilization of the metastable state and res- 
onant domain wall propagation, that are not observed 
in equilibrium, but are likely to characterize a broad 
class of actual systems dominated by dynamic impuri- 



ties. Surprisingly, these far-from-equilibrium phenomena 
can be understood using nucleation theory based on a 
nonequilibrium potential at a mesoscopic/cluster level. 
This is possible because the relevant excitations driving 
the metastable-stable transition result from a competi- 
tion between a surface and a bulk. This suggests that our 
approach may prove useful for studying metastability in 
many other nonequilibrium systems where a surface-bulk 
competition plays a significant role. 
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